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SUMMARY 

In  order  to  use  the  visual  evoked  potential  as  a  real-time  monitor  of  the  state  of  the  visual  system,  a 
monitoring  technique  was  developed  to  maximize  the  signal-to-noise  (S/N)  ratio.  System  response 
parameters  (co-power,  phase  angle,  and  coherence)  are  estimated  for  each  subject  from  an  ensemble 
of  FFT's  during  a  pre-testing  period.  These  parameters  are  then  used  to  calculate  a  weighting  func¬ 
tion  which  is  used  in  real-time  to  perform  amplitude  normalization,  coordinate  rotation,  and  op¬ 
timal  weighting  of  the  terms  of  individual  FFT's.  These  terms,  when  summed,  produce  a  simple 
real  variable  with  maximal  S/N  and  an  expected  value  of  100  for  normal  vision  and  zero  for  no 
vision  or  blackout.  This  technique  has  been  implemented  with  an  FFT  signal  analyzer  and  a  desk¬ 
top  computer.  Experimental  results  for  six  subjects  indicate  that  a  useable  measure  may  be 
achieved.  While  improvements  in  the  method  are  required  if  we  are  to  have  an  effective  real-time 
monitor  for  visual  functioning,  such  improvements  are  not  only  possible,  but  feasible.  Optimization 
of  variables  such  as  lead  placement  and  stimulus  timing  promise  the  required  improvement. 


INTRODUCTION 

Modern  high-performance  aircraft  can  produce  and  sustain  levels  of  acceleration  {+Gz)  greater  than 
the  pilot's  physiologic  tolerance  limits.  These  limits,  as  a  joint  function  of  -i-Gz  parameters  and 
numerous  protective  devices  and  methods,  are  studied  on  human  centrifuges  such  as  the  Dynamic 
Flight  Simulator  at  the  Naval  Air  Development  Center. 

Subjects  can  experience  loss  of  consciousness  (LOC)  during  exposure  to  sufficiently  high  -t-Gz 
acceleration.  LOC  is  usually  preceeded  by  visual  symptoms:  decreased  visual  sensitivity,  dimming 
of  the  visual  field,  peripheral  light  loss  (PLL),  and  central  light  loss,  ending  in  complete  light  loss  or 
blackout.  PLL  is  commonly  used  as  an  end  point  indicator  for  centrifuge  runs  [1] . 

Real-time  monitoring  of  PLL  requires  an  active  response  by  the  subject.  In  contrast,  the  visual 
evoked  potential  (VEP)  [2] ,  which  also  reflects  the  integrity  of  the  visual  system,  requires  only  a 
relatively  passive  task  of  viewing  a  visual  stimulus. 

Real-time  monitoring  of  VEP  as  a  measure  of  visual  functioning  on  the  human  centrifuge  must  meet 
certain  criteria,  and  stay  within  certain  constraints.  The  most  important  are  that  the  response-time 
to  loss  of  vision  must  be  less  than  approximately  four  seconds,  with  a  low  false-alarm  rate  and  a 
very  low  probability  of  not  detecting  (and  acting  upon)  a  complete  loss  of  vision.  The  method  must 
not  require  extended  subject  preparation  time,  e.g.,  a  many-electrode  montage.  The  required  com¬ 
puter  capacity  must  be  less  than  a  full  mini-computer  and  the  non-real-time  computation  must  not 
significantly  extend  the  normal  1-2  minute  rest  periods  between  centrifuge  runs. 

The  available  real-time  methods  (e.g.,  [3,  4,  5] )  did  not  seem  to  be  directly  applicable  when  con¬ 
sidered  in  terms  of  these  criteria  and  constraints.  However,  such  methods  did  provide  the  broad 
strategy  for  attempting  a  solution. 

To  provide  a  continuous  flow  of  data,  we  chose  to  utilize  the  steady-state  VEP  [3] .  A  monitoring 
method,  staying  within  these  computer  capacity  and  time  constraints,  was  developed  [6] .  This 
method  was  implemented  on  a  desk-top  computer  (HP  9825T)  and  an  FFT-based  low-frequency 
dual  channel  signal  analyzer  (HP  3582A),  and  tested  with  six  subjects  under  static  (non-centrifuge) 
conditions.  The  results  proved  encouraging  and  potentially  useful  in  general  VEP  studies  [7] . 
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METHOD 


OVERVIEW 

The  method  requires  an  initial  parameter  estimation  run,  relating  system  output  (VEP  &  EEC)  to 
input  (stimulus).  Based  on  an  ensemble  of  FFT's,  calculation  provides  co-power,  cross-power,  phase 
angle,  and  coherence  [8] .  These  are  used  to  calculate  a  weighting  function  to  be  applied  to  single 
FFT's  in  real-time,  to  produce  a  single  variable  which  may  be  labeled  "Percent  Vision".  The  weight¬ 
ing  function  performs  magnitude  normalization,  rotation  of  axes  such  that  the  VEP  components 
become  positive-real,  and  optimal  weighting.  When  summed  together,  these  optimally  weighted 
components  become  the  variable  "Percent  Vision".  Under  the  restrictive  assumptions  of  our  model, 
this  variable  should  have  a  maximal  S/N,  a  predicted  standard  deviation,  an  expected  value  of  100 
under  conditions  of  normal  vision,  and  0  under  conditions  of  blackout. 


PARAMETER  DETERMINATION 

In  order  to  determine  the  basic  parameters  of  the  process  we  collected  simultaneous  ensembles  of 
input  (stimulus)  and  output  (response)  FFT's  from  the  dual  channel  frequency  analyzer.  Here  we 
designate  the  input  as  X(f)  and  the  output  as  Y(f)  where  Y(f)  is  actually  VEP-t-EEG.  The  computa¬ 
tion  of  the  other  parameters  is  as  follows:  Input  auto-power  spectrum 

Gxx  (f)  =  X(f)  •  X*  (f)  ( 

output  auto-power  spectrum 

Gyy  (f)  =  Y(f)  •  Y*  (f)  ( 

and  cross-power  spectrum 

(f)  =  Y(f)  •  X*  (f)  ( 

Then,  using  the  ensemble  averages,  the  coherence  function 


Gxx  (f)  •  G,,  (f) 

SIMPLIFIED  PARAMETER  DETERMINATION 

We  are  concerned  only  in  regularities  in  the  data,  not  with  veridicality.  In  real-time,  vve  want  to 
process  a  raw  FFT  to  produce  a  single  optimized  variable.  At  the  level  of  the  FFT,  this  means  that 
we  do  not  want  corrections  for  phase-shift  and  pass-band  ripple  produced  by  the  anti-aliasing  and 
decimation  filters  in  the  signal  analyzer.  Similarly  we  used  a  uniform  window  function  for  the  time 
function  measurement. 

By  careful  design  and  control  of  the  stimulus,  a  number  of  advantages  and  simplifications  occur. 

For  the  VEP,  using  a  strobe  light  driven  by  a  pulse  train  or  trains,  each  response,  whether  funda¬ 
mental  or  harmonic,  is  an  integral  multiple  of  the  stimulus  fundamental,  which  is  chosen  to  be  an 
integer  multiple  of  the  bin  spacing.  This  eliminates  ''picket-fence"  effects  and  places  the  response 
at  specific  harmonically  related  spectral  points. 
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The  stimulus  and  measurement  processes  are  slaved  such  that  each  stimulus  time-frame  is  identical 

[Xj  (t)]j  =  [Xj  (t)] 


% 


% 


* 


# 


Since  the  stimulus  is  thus  constant,  then  for  our  relativistic  approach  in  the  detection  process  we 
define  the  input  amplitude  spectrum  as  uniform  and  unity  magnitude  with  zero  phase  angle  over 
the  band  of  interest.  This  eliminates  the  need  to  calculate  the  FFT  for  the  input  and  greatly  sim¬ 
plifies  the  calculations  required  later  on.  Thus 


X(f)  =  1  (5) 

and  therefore 

G,,  (f)  =  1.  (6) 

Then,  by  substituting  into  equation  (3): 

Gy,  (f)  =  Y(f)  (7) 

and,  finally 

72  =  |Y(f)|2/Gyy  (f).  (8) 


Also,  note  that  Y(f)  directly  specifies  the  phase  angle  of  the  response. 

Since  a  response  will  be  obtained  at  the  fundamental  and  possibly  harmonic  frequencies  of  the 
stimulus  X(f),  the  FFT  need  only  be  processed  at  those  frequencies.  Thus,  in  our  application,  with 
a  0.8  Hz  spacing  of  a  128  point  FFT  over  the  0-100  Hz  band,  using  a  mixed  pulse  train  of  stimuli 
(12.8  and  16  Hz)  synchronized  with  a  missing  fundamental  of  f^  =  3.2  Hz,  we  analyzed  at  every 
fourth  non-zero  point,  for  a  total  of  31  points. 

To  summarize,  our  simplified  computational  requirements  consist  of: 

1 )  For  each  time  frame  of  the  ensemble,  compute  Y(f)  and  decimate  to  include  only  multiples 
of  the  missing  fundamental:  YlNf^)  where  N=1,...,  31.  For  each  time-frame,  compute  auto¬ 
power  spectrum,  G^y  (Nf^). 

2)  Compute  ensemble  averages,  YINf^)  and  Gyy  (Nf^). 

3)  Compute  coherence. 


7^  (NfJ  =  |Y(NfJ|2/Gyy  (NfJ  (9) 

Unless  the  number  of  frames  in  the  ensemble  (N^)  is  quite  large,  coherence  should  be  corrected  for 
bias.  Thus,  setting  l/N^  as  the  threshold  criterion. 

If  72  (Nfjj )  <  1/Nf,  then  72  =  0  and 

if  7^  (Nf„)>1/Nf,then 

72  -  1/Nf 
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(All  further  reference  to  coherence  assumes  this  correction). 

CREATION  OF  NORMALIZED  REAL  VARIABLES 

The  corrected  coherence  permits  us  to  partition  the  power  (variance  about  the  origin)  into  signal, 

S=72Gyy  (11) 

and  noise. 


N  =  (1  -7'')Gyy  =  Gyy  -S  (12) 

components.  Since  the  noise,  which  is  predominantly  background  EEG,  is  phaseless,  its  variance 
(which  is  produced  by  squaring  and  summing  two  zero  mean,  equal  variance,  orthogonal,  linear 
variables)  is  equal  to  N/2  along  any  arbitrary  axis.  That  is,  phase  sensitive  detection  rejects  Va  of 
the  incoherent  noise.  Th^n,  defining  a  new  axis  system  where  the  positive  real  axis  is  aligned  with 
and  of  the  same  sense  as  Y(NfQ),  the  new  real  component  can  be  determined  by  defining  a  rota¬ 
tional  weight,. 


Y(NfJ 
W,  =— - 

|Y(NfJ| 

and  calculating  the  function, 

(ReWi)  (Re  Y)  -i-  (Im  W^)  (Im  Y). 

Using  the  expected  standard  deviation  about  the  origin  we  derive  a  normalizing  weight, 
W2  =  l/(  (72  +’72(1-72))  Gyy)’/>  =  l/(y2(l+72)Gyy)’''' 

Then,  these  weights  may  be  combined  as. 


(13) 


(14) 


(15) 


W3  =  Re(WiW2),  lm(WiW2).  (16) 

We  may  then  calculate  real,  positive  mean,  unity  power  variables,  using  single  FFT's  in  real-time, 

R  =  (ReW3)  (ReY)  +  (lmW3)  (ImY)  (17) 

which  have  the  enhanced  coherence. 


7i  = 


7^  +  ’72(1-72) 


272 

1+7-2 


(18) 


(All  further  use  of  coherence  assumes  7I ). 

By  combining  decimation  across  frequency  with  phase  sensitive  detection,  we  have  improved  overall 
S/N  by  a  factor  of  eight.  By  conversion  to  normalized,  positive  mean,  real  variables,  we  are  in  a  posi¬ 
tion  to  optimally  weight  and  sum  these  variables  for  further  S/N  improvement. 
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OPTIMAL  WEIGHTING 


In  the  most  general  form,  the  variance  of  a  weighted  sum  of  variables  involves  a  full  variance- 
covariance  matrix  [9] .  It  was  decided  that,  while  multiple  regression  [10]  and/or  linear  discrim¬ 
inant  [11]  models  could  be  applied,  the  required  computation  would  exceed  our  capacity-time 
constraints.  If  we  assume  that  all  of  the  structure  of  the  covariance  of  the  variables  is  due  to  unity 
correlation  of  signal  components  with  each  other  and  zero  correlation  of  noise  components  with 
each  other  and  with  signal,  then  a  step-wise  process  may  be  used.  Here,  two  variables  are  combined, 
the  compound  is  then  combined  with  the  third,  the  resulting  compound  is  then  combined  with  the 
fourth,  and  so  on.  Computational  requirements  increase  as  a  linear  function  of  the  number  of 
variables  to  be  combined. 

For  two  variables  [9] ,  setting  W,  =1 ,  we  may  solve  for  relative  weight  (W,  =1,  W2=W) : 

■H2ri2Waia2.  (19) 


The  variables  which  we  are  weighting  and  summing  are  real  and  each  have  two  components  with 
variances  (powers)  Sp,  and  Np, ,  such  that  Sp,  +  Np,  =1.  The  signal  variance  of  the  resultant  com¬ 
pound  is, 

Src  =  Sri  +  W2Sr2 +2W(SriSr2)’''’  '  (20) 


the  standard  deviation  is. 


and  the  noise  variance  is. 


(21) 


Nrc  =  Nri  +  W2Nr2  =  (1-Sri)  +W2  (1-Sr2).  (22) 

To  find  the  weight  W  producing  the  maximum  Sr^/Nr^  ,  we  differentiate  it  with  respect  to  W  and 
set  the  result  equal  to  zero. 


0  =  [(Sr2-1)  (Sr iSr2 )’''’]  W2 

+  [Sr2-Sri]W  (23) 

+  [(1-Sr  1  )  (Sr  ,  Sr2  )’'^’  ]  • 

This  quadratic  equation  is  then  solved  for  W. 

Some  rather  tedious  change-of-scale  arithmetic  is  involved  in  getting  from  W  as  delivered  by  equa¬ 
tion  (23)  to  the  actual  weighting  function  to  be  used  for  optimization.  In  our  application,  we  take 
the  first  of  the  compressed  spectral  points  and  set  W,  =0  if  7^  <0.15,  and  arbitrarily  set  W,  =7^  if 
not.  This  defines  values  "Src ''=W^7^ ;  "Npg ''=W^  (1-7^;  and ''7^''=7^.  The7pg  is  passed  as 
Sr,  to  the  optimization  equation  along  with  72  (if  >  0.1 0)  as  Sr2-  Th 
which  assumes  that  S,  +  N,  =  1,  is  rescaled  to. 


ie  delivered  weight,  , 


^2  (2rci  ■'■^rci)^ 


(24) 
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Wj  is  then  used  to  produce  the  first  genuine  compound.  Since  signals  add  as  standard  deviations 
(eq.  21 )  and  noises  as  variances  (eq.  22), 

2rc2  “  ^^RCi 

^RC2  ~  ^RCi  W^(1-S2),  and 
T'rC2  ~  ^RC2/^^RC2  NrC2^- 

This  step-wise  process  continues  until  all  useable  points  are  exhausted. 

Final  values  of  Src  and  Nr^  are  used  to  rescale  the  entire  weighting  function,  to  produce  an 
expected  value  of  100%.  This  weighting  function  is  combined  with  the  normalization  and  rotation 
weights  to  produce  the  weighting  function  for  use  in  real-time  runs.  The  final  value  of  I'Trc 
used  to  set  confidence  limits  for  the  variable  "Percent  Vision",  which  is  computed  and  displayed 
in  real-time. 

The  final  7^^.  also  provides  a  simple  figure-of-merit  for  evaluation  of  any  change  in  experimental 
method.  Changes  in  electrode  placement,  changes  in  stimulus  parameters  (e.g.,  intensity,  frequency 
composition,  spatial  structure),  or  any  other  changes  in  method  may  be  evaluated. 

The  method  is  presently  being  used  to  process  only  a  single  channel  of  information:  the  VEP  mea¬ 
sured  at  the  mid-line  (0^).  Separate  processing  for  right  and  left  hemispheres  (0^ ,  Oj)  should  pro¬ 
vide  a  significant  improvement,  and  still  stay  within  reasonable  bounds  for  required  computing 
capacity/time.  The  algorithm  would  remain  the  same,  treating  the  two  channels  as  a  single  set. 

EXAMPLE 

Figure  1  shows  a  specimen  real-time  run  for  a  subject  with  low  predicted  and  obtained  standard 
deviation  [7] .  Here,  the  stimulus  was  occluded  for  three  time-frames  in  the  latter  part  of  the  run. 
The  dashed  line  in  the  figure  shows  the  result  of  a  two-point  smoothing  which  provides  a  margin  of 
protection  against  a  false  negative  determination  at  the  expense  of  an  additional  1.6  second  delay 
of  the  decision. 
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Figure  1.  Specimen  Real-Time  Run  (from  [7]). 


7 


NADC-82206-60 


REFERENCES 

[1]  M.  M.  Cohen,  "On  the  additivity  of  acceleration  protection",  in  Aerospace  Medical  Assoc. 
Annual  Meeting  Preprints,  May  1981,  pp.  9-10. 

[2]  D.  Regan,  Evoked  Potentials  in  Psychology,  Sensory  Physiology  and  Clinical  Medicine, 
London,  Chapman  and  Hall,  1972. 

[3]  D.  Regan,  "Steady-state  evoked  potentials",  J.  Opt.  Soc.  Amer.,  vol.  67,  1977,  pp.  1475- 
1489. 

[4]  T.  F.  Collura,  "Real-time  evoked  potential  correlates  of  directed  attention",  in  Proc. 
ACEMB,  Oct.  1978,  p.  10. 

[5]  R.  W.  Sencaj,  J.  I.  Aunon,  and  C.  D.  McGillem,  "Discrimination  among  visual  stimuli  by 
classification  of  their  single  evoked  potentials",  Med.  Biol.  Eng.  Comput.,  vol.  17,  1979, 
pp.  391-396. 

[6]  J.  G.  Nelson,  J.  J.  Palumbo,  and  M.  M.  Cohen,  "Use  of  the  coherence  function  in  analysis 
of  steady-state  visual  evoked  potential  (VEP)",  in  Research  Program  Review,  Aircrew 
Physiology,  Naval  Air  Systems  Command  Research  and  Technology  Group,  Oct.  1980, 
pp.  15-24. 

[7]  J.  G.  Nelson  and  L.  Hrebien,  "Real-time  monitoring  of  the  visual  evoked  potential",  in 
Aerospace  Medical  Assoc.  Annual  Meeting  Preprints,  May  1982,  pp.  169-170. 

[8]  J.  S.  Bendat  and  A.  G.  Piersol,  Random  Data:  Analysis  and  Measurement  Procedures,  New 
York,  Wiley-Interscience,  1971. 

[9]  J.  P.  Guilford  and  B.  Fruchter,  Fundamental  Statistics  in  Psychology  and  Education  (6th 
Ed.),  New  York,  McGraw-Hill,  1977. 

[10]  N.  P.  Draper  and  H.  Smith,  Applied  Regression  Analysis  (2nd  Ed.),  New  York,  Wiley,  1981. 

[11]  D.  F.  Morrison,  Multivariate  Statistical  Methods  (2nd  Ed.),  New  York,  McGraw-Hill,  1976. 


« 


REPORT  NO.  NADC-82006-60 


DISTRIBUTION  LIST  (continued) 

No.  of  Copies 

Dr.  Jerry  I.  Nelson,  Assistant  Professor  of  Experimental 

Opthalmology,  N.Y.  University .  1 

Royal  Air  Force  Institute  of  Aviation  Medicine, 

Farnborough,  U.K .  1 

Defense  &  Civil  Institute  of  Environmental  Medicine, 

Downsview,  Canada .  1 

German  Air  Force  Institute  of  Aerospace  Medicine, 

F'urstenfeldbruck,  West  Germany .  1 

Institute  of  Aviation  Medicine,  Oslo,  Norway .  1 


REPORT  NO.  NADC-82006-60 


DISTRIBUTION  LIST 


No.  of  Copies 


Director,  Defense  Technical  Information  Center . 

Commander,  Naval  Air  Systems  Command . 

(4  for  NAVAIRDEVCEN  Liaison  Office) 

(3  for  AIR-310C) 

(2  for  AIR-340B) 

(1  for  AIR-531B) 

Chief  of  Naval  Material . 

(1  for  MAT-05) 

Commanding  Officer,  Naval  Medical  Research  &  Development  Command.. 
(1  for  NMRDC-44) 

Chief,  Bureau  of  Medical  &  Surgery . 

(1  for  NM&S  3C1) 

Chief  of  Naval  Operations . 

(1  for  N0P-05H) 

(1  for  N0P-09E) 

Chief  of  Naval  Research . 

(1  for  ONR-440) 

(1  for  ONR-441) 

(1  for  0NR-441NP) 

(1  for  ONR-442) 

Commander,  Naval  Safety  Center . 

Commanding  Officer,  Naval  Aerospace  Medical  Research  Laboratory... 

Superintendent,  Naval  Postgraduate  School . 

Commanding  Officer,  Naval  Health  Research  Center . 

Personnel  Research  &  Development  Center. 

Center . 

Biodynamics  Laboratory . 

Submarine  Medical  Research  Laboratory... 

Training  Equipment  Center . 

Scientific  Research  (AFSC)/NL . 

Air  Force  Aerospace  Medical  Research  Laboratory . 

U.S.  Air  Force  School  of  Aerospace  Medicine . 

U.S.  Army  Aeromedical  Research  Laboratory . 

FAA  Civil  Aeromedical  Institute . 

NASA  Ames  Research  Center . 

NASA  Johnson  Space  Center . 

Dr.  Lenard  Matin,  Visual  Science  Laboratory,  Columbia  University. . 
ron  Lewis,  School  of  Applied  &  Engineering  Physics, 

rnell  University . 

A.  TenjjiiflPlfciBahill ,  Biomedical  Engineering  Program, 

Carne^fe  Mel#on  Universi^.  —  . . 

.  Bruce  A.  E*enstein,  lUik  Department*! ectrical 
Engineerinifand  Dr.  B^  (liaral,  D|^e\  Uni versi 
Dr.  Mike  L.  Jalbarsht»jFofe|bor  o^^tha  jiology  and, 

Enginea^ng,  Duke, 


12 

10 


Commanding 
Commander, 
Command i ng 
Command i ng 
Commanding 
Air  Force 


Officer,  Naval 
Naval  Air  Test 
Officer,  Naval 
Officer, 
Officer, 

Office  of 


Naval 

Naval 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

2 

1 

1 

1 

2 

1 

1 


il  &  Computer 

;ityj# . ^ 

md^omedical  M 

/  ^ 


